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It has proven a significant challenge to experiment and phenomenology to extract precise 
values of the nucleon sigma terms. This difficulty opens the window for lattice QCD 
simulations to lead the field in resolving this aspect of nucleon structure. Here we report 
on recent advances in the extraction of nucleon sigma terms in lattice QCD. In particular, 
the strangeness component is now being resolved to a precision that far surpasses best 
^-i ■ phenomenological estimates. 

i> : 

in . 1. NUCLEON SIGMA TERMS 
l> ; 

The nucleon sigma terms are important quantities in resolving the dynamics of QCD, 
where they help to understand the role of explicit chiral symmetry breaking in the mass 
Q\ [ of the nucleon. Indeed, the strangeness component plays a unique role in that is purely 
a vacuum polarization effect, analogous to the Lamb shift in QED. Further, the strange 
quark is light enough that it probes the nonperturbative distance scales of QCD. Beyond 
^ . the fundamental importance to nucleon structure, the strange quark condensate in the 
nucleon is of significant interest in studies of the QCD phase structure at large baryon 
density [1] and in constraining predicted cross sections for dark matter detection [2113] . 

The sigma terms of the nucleon are defined by the scalar form factors in the limit of 
vanishing momentum transfer, a q = m q (N\qq\N), where q denotes the quark flavour of 
interest. The light-quark sigma term (or pion-nucleon sigma term), <Ji = me(uu + dd) 
(with mi = (^u+ m ci)/2), is related to n-N scattering through a chiral low-energy relation 

= <j£ = £.Jv — Ajt — A a , (1) 

where the experimental input required is the Born-subtracted, isoscalar ttN scattering 
amplitude evaluated at the (unphysical) Cheng-Dashen point, The smallness of the 

up and down quark masses ensures that both the remainder term, A# [5T7] , and the form 
factor correction, A a [8], can be reliably estimated. An early analysis of experimental 
results [H] led to the value = 45 ± 8 MeV |S] , while a more recent analysis suggested a 
higher value, 64 ± 7 MeV [TO]. 

Unlike the light quarks, the strange quark is too heavy to invoke a low-energy relation. 
Therefore a s has generally been estimated by studying the breaking of SU(3) within the 
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baryon octet [Tfll|ll2j . Here, the observed spectrum has been used to derive a constraint 
on the non-singlet combination 

(T = m e (N\uu + dd- 2ss\N) , (2) 

where chiral effective field theory studies lead to a value cr = 36 ± 7 MeV [12] , building 
on earlier estimates by Gasser [IT]. The difference between ctq and the extracted <jg 
then defines the strange-quark contribution, <jt — o"o = 2^ L a s . Equivalently, the relative 
strangeness component is often discussed in terms of the parameter y, 

q _ 2(N\aa\N) _ 

a e (N\uu + dd\N) V ' 1 ' 

Taking the above pion-nucleon sigma term values and the estimated <To gives values y ~ 
0.2 ± 0.2 and 0.44 ± 0.13, and using a strange-to-light quark mass ratio of ~ 25 gives 
strangeness sigma terms, a s = 110 ± 130 and 350 ± 120 MeV, respectively. This approach 
clearly leads to a result for a s which is very sensitive to the precise value of a£. We also 
note that even with perfect n-N data to better constrain <jg, the 7 MeV uncertainty in ao 
alone leads to a ~ 90 MeV uncertainty in a s . This limitation clearly opens the way for 
lattice QCD to offer significant improvement. 

2. LATTICE CALCULATIONS 
2.1. Light quarks 

In early lattice simulations, limited computing power forced the use of heavy quarks 
and the quenching of the light quarks in the theory [T3~yT4"] . The first direct calculations 
of the disconnected contributions (though still on a quenched gauge ensemble), including 
an estimate of the strange quark contributions, were performed in the work of Fukugita et 
al. [15] and Dong et al. [16]. In these works, the 7r-N sigma was found to be consistent 
with observation — see Figure [JJ for a summary of results. 

The first determination using (2 flavours of) dynamical quarks found a value for ag 
that is smaller than phenomenological expectations [17] . A notable feature of this cal- 
culation is that consistent results were demonstrated between the direct calculation of 
the disconnected loops and that extracted from the fit to the nucleon mass, using the 
Feynman-Hellman relation^ 

dM N 

a « = m ^ ■ (4) 

While this calculation appeared to yield a value that was rather low, it was soon identified 
that the sigma term is particularly sensitive to chiral extrapolation [T9] . After this time, 
most extractions of a£ have come by application of the Feynman-Hellman relation to the 

It was common in the early calculations that the quoted error only reflected the sta- 
tistical uncertainty in the underlying extraction. The work of Procura et al. [22] made 
a detailed investigation of the systematic errors associated with the input parameters of 



1 Commonly the quark masses in this relation are replaced by the meson mass squared using the Gell- 
Mann-Oakes-Renner relation. 
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Figure 1. Summary of lattice QCD calculations of the light-quark sigma term (left panel) 
and strangeness sigma term (right panel). Triangles, squares and circles denote quenched, 
2-flavour and 2+1-flavour lattices, respectively. We caution that different calculations 
have invested differing degrees of effort into quantifying their error estimate. 



their functional form. Contrasting this with their earlier work [21] , where only the statis- 
tical error was quoted, is a useful guide as to the potential size of systematic effects. Of 
course, various systematics will continue to be better controlled as the simulation results 
improve. 

While it is not obvious from Figure [H(a) that the lattice results are showing convergence 
- without taking bias — there does appear to be consensus that the lattice results are 
consistent with the two phenomenological extractions discussed above. 

2.2. Strange quark 

Early lattice estimates of the strangeness sigma term supported the rather large values 
that had been inferred from phenomenology [ToTTE] . Going beyond the quenched ap- 
proximation, the 2-flavour dynamical calculations of the SESAM Collaboration [17] also 
supported the larger strangeness values, though with substantial uncertainty. Another 
early calculation highlighted some of the challenges of resolving these matrix elements 
and concluded that the strangeness component is most likely rather small |18j . 

While the strange quark in the 2-flavour simulations is essentially quenched, it is inter- 
esting to consider the light-quark disconnected diagrams, which are truly dynamical. One 
can consider the ratio of the disconnected to connected (or sea to valence) contributions 
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to the sigma term, 

(N\uu + dd\N) 

disc /-\ 

d/c = (N\uu + dd\N) con 1 j 

In the quenched simulations a relatively large value for this ratio was observed, with 
R d / C = 2.23(52) [15] and 1.79(7) [IB]. We note that the effect of unqenching appears to 
reduce this ratio, with the corrsponding ratio found to be Rd/c — 1-26(57) in the dynamical 
simulation [IT]. This indicates a decrease in the relative strength of the disconnected 
component. 

Extending this discussion, a modern determination of this ratio by JLQCD [24] has 
found Rd/c = 0.41(5) (taking PQ-b as the "best estimate") — suggesting even further re- 
duction of this ratio as results better describe the physical regime. This study constrained 
the nucleon mass using light-quark masses taking values over the range m s /6-m s (or pion 
masses ~ 290-750 MeV). 

In a SU(3) symmetric (m u = = m s ) 3-flavour dynamical calculation, the discon- 
nected contributions of each of the three quarks are identical, satisfying the identity 

2(N\ss\N) disc SU ^ 3) (N\uu + dd\N) disc . (6) 

In a 2-flavour simulation, the literal ss contribution is zero, yet it could be reasonable 
to adopt the above relation to infer the ss from a 2-flavour calculation. Indeed this is 
the approach taken by [24J, where at the SU(3) symmetric point one can infer y ~ 0.09. 
Further taking the 2-flavour result for the increase in both the connected and disconnected 
light-quark terms as the light-quarks are reduced to the physical point reduces the estimate 
of y to ~ 0.03 [2l] (assuming that the strangeness term is only weakly dependent on the 
light-quark masses). Indeed, when combined with their extraction for the o>, this leads 
to a significantly reduced value of the strangeness sigma term <j s = 20 ± 12 MeV. 

The major cause of this reduction compared with the previous estimates comes from 
the observation that the derivative of the nucleon mass ( is quite strongly dependent 

v UTYlsea 

on the sea-quark mass (m sea ). Indeed, this derivative was observed to be enhanced by 
a factor of ~10 in going from the strange-quark mass to the light-quark mass. Beyond 
the numerical improvements in the underlying lattice simulations, this study has also 
incorporated the consequences of chiral symmetry breaking, where there are known to be 
substantial chiral logarithms in these derivatives — leading to the observed enhancement 
in the light-quark domain. 

These authors have recently reported results on testing their method against a direct 
calculation of the disconnected strangeness component, finding compatible results |27j . 
Further, the analysis of the JLQCD Collaboration is being extended to 2+1 flavours of 
dynamical quarks, where preliminary results also support their 2-flavour findings |28j . 

Another opportunity to study the strange-quark mass dependence came with a series of 
calculations of the baryon spectrum in 2+1-flavour dynamical simulations which emerged 
in 2008 [29 30 3 11132] (after an early venture in 2001 [33J). These new 3-flavour results 
allowed tests of the SU(3) chiral expansion at the lattice quark masses [29.25 26J. Indeed 
the poor convergence of the SU(3) expansion, already recognised at the physical quark 
masses [T2~|ITT] . was rediscovered on the lattice. These three chiral-lattice papers took 
three different approaches to deal with this problem. PACS-CS [26] simply abondoned 
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the SU(3) formulation; LHPC |29j sacrificed the coefficients of the chiral logarithms in 
order to stabilise the fits; we [25J introduced a single new fit parameter that acts to 
separate the low- and high-energy contributions to the chiral loop integrals. 

The results of Ref. [25] built on the established benefits of using finite-range regulari- 
sation, which has previously been shown to dramatically improve the SU(3) convergence 
[TT1I34"] and offer robust chiral extrapolation for lattice QCD [20.35J. The analysis of 
Ref. [25J demonstrated the ability to extrapolate in both the light- and strange-quark 
masses to accurately reproduce the physical octet baryon spectrum, and also extrapolate 
to the heavier quark masses (not used in the fits) as calculated on the lattice. The results 
were found to be consistent with two different lattice discretizations, which also allowed 
an estimate of potential discretization effects. 

Using the Feynman-Hellman relation, we found cr^ to be compatible with phenomenolog- 
ical and lattice estimates. And the strangeness sigma term was identified to be relatively 
small, supporting the 2-flavour calculation of JLQCD [24"] . 

Another recent 2+1-flavour dynamical result for a s has been reported by Toussaint 
and Freeman [36], using an application of the Feynman-Hellman relation on the nucleon 
correlator. This technique substantially differs from the approaches based on fitting the 
the nucleon mass and therefore provides a strong independent test. The result of Toussaint 
and Freeman [SB] is marginally higher than the values reported by JLQCD [23] and by 
us [25] . Nevertheless, as seen in Figure W[b), the modern lattice results for cr s are in 
agreement that the size of the strangeness sigma term is substantially smaller than has 
been previously thought. 

3. IMPACT 

As described above, the dominant uncertainty in knowledge of the scalar quark cou- 
plings has lies in the strangeness component. The current generation of lattice QCD 
calculations has primarily resolved this limitation by providing stringent new limits on 
the strange quark sigma term. As a result of this improvement, the predicted cross sec- 
tions for models of dark matter can now be revised. In a recent paper, Giedt et al. [37] 
have investigated the impact of the new lattice QCD results on the cross sections for a 
class of super symmetric models. The poorly known strangeness cx-term had previously 
resulted in cross sections which varied by roughly an order of magnitude [3] . The uncer- 
tainties have been dramatically reduced by the lattice calculations and importantly, one 
now has significant discrimination power between the ensemble of models considered. 

This work was supported by the U.S. Department of Energy, Office of Nuclear Physics, 
under Contract No. DE-AC02-06CH1I357; and the Australian Research Council. 
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